Anisotropy of the upper critical field in superconductors with anisotropic gaps. 
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The upper critical field Hc2 is evaluated for weakly-coupled two-band anisotropic superconductors. 
By modeling the actual bands and the gap distribution of MgB2 by two Fermi surface spheroids 
with average parameters of the real material, we show that Hc2.ab/ Hc2,c increases with decreasing 
temperature in agreement with available data. 

PACS numbers: 74.60.Ec, 74.20.-z, 74.70.Ad 



The anisotropic Ginzburg-Landau (GL) equations, 
derived for clean superconductors with arbitrary gap 
anisotropy by Gor'kov and Melik-Barkhudarov led to 
a common practice of characterizing materials by a single 
anisotropy parameter defined as £,a/£.c = ^c/^a {£, is the 
coherence length, A is the penetration depth, and a, c are 
principal directions of a uniaxial crystal of the interest 
here). Formally, this came out because in the GL do- 
main, the same "mass tensor" determines the anisotropy 
of both ^ (of the upper critical fields i?c2) and of A. 

At arbitrary temperatures, however, the ratios of Hc2S 
and of A's are not necessarily the same. We demonstrate 
below that in materials with anisotropic Fermi surfaces 
and anisotropic gaps, not only i?c2,a/-ffc2.c may strongly 
depend on T, but this ratio might differ considerably 
from Xc/Xa at low T's. Our arguments are based on 
the weak-coupling model of superconductivity for simple 
Fermi surfaces and gap anisotropics; as such they are 
at best qualitative. Still, being applied to MgB2, they 
provide satisfactory description of existing data for the 
Hc2 anisotropy. 

There are two different approaches in literature in de- 
scribing macroscopics of the "two-gap" superconductiv- 
ity of MgB2. One of these treats the anisotropy of in- 
teraction by introducing a coupling matrix pij for intra- 
and interband pair transfer. Various relations between 
the matrix elements yield variety of macroscopic conse- 
quences. One of the realizations of the model considers 
two order parameters with two distinct phases which may 
lead to various static ^ and dynamic Q effects. Cer- 
tain relations between elements pij provide the experi- 
mentally observed gaps and their ratio. The approach 
adopted in this paper considers the gap on two Fermi 
sheets just as a particular case of the gap anisotropy, 
the gap ratio being the experimental input parameter. 
Within this scheme, there is only one complex order pa- 
rameter ^P, a single critical temperature Tc is built in, and 
the number of input parameters needed for calculations 
is small. The resulting Hc2^s of these two approaches are 
similar ||^ . Our choice is dictated by theoretical simplic- 



ity, rather than by experimental necessity. 

Below, the general scheme of calculating Hc2(T) 
based on the Eilenberger method Q is outlined for 
both the gap and the Fermi surface being anisotropic 
0. Then, we estimate the Hc2 anisotropy of MgB2 
by modeling the four-sheet Fermi surface as calculated 
in Refs. ^, |^, |l^ with the gap distribution having 
two sharp maxima, by two distinct Fermi sheets F12 
with gaps Ai^2, each being constant within its sheet. 
Since the actual band structure enters macroscopic 
superconducting parameters via Fermi-surface averages, 
we further model the sheets -^1^2 by two spheroids with 
average Fermi velocities close to the band-structure 
generated values. As a result we obtain qualitatively 
(and - given the spread of existing data - quantitatively) 
correct behaviors of Hc2{T) for both principal directions. 

A clean weakly coupled superconductor at Hc2(T) is 
described within the quasiclassical Eilenberger scheme by 
the anomalous Green's function /(tj, v, r), which satisfies 
an equation: 



{2uj + vU)f = 2Alh. 



(1) 



Here, v is Fermi velocity; lo = TTT{2n + are the 
Matsubara frequencies; 11 = V -I- {2TTi/(f>Q)A with the 
vector potential A and the flux quantum 0o- The Eilen- 
berger function g describing normal excitations, is unity 
at Hc2 (T) . Equation (|l]) holds for any Fermi surface and 
gap anisotropics. 

The gap function A satisfies the self-consistency equa- 
tion: 



A(r, v) = 27rTiV(0) V {V{v, v' )f{uj, v' , r))„, 



CJ>0 



(2) 



Here, (...) denotes the average over the Fermi surface. 
The pair scattering potential V is assumed factorizable, 
V{v,v') — Vq il,{v) il,{v' ) , with the function n{v) nor- 
malized so that 
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(3) 
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Further, we look for A(r,T;t!) = ^{r,T)n{v) ||Tl[, the 
form implying a one-component order parameter. Then, 
after excluding Vq with the help of the BCS formula for 
Tc, one arrives to Eq. (0) of the form: 



27rT T 



°° ,T. 



w>0 



V huj 



(4) 



(for more details and references see, e.g., Ref. p^ . 
The linear Eq. (|^) is inverted: 
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(5) 



In the GL domain, the gradients 11 ~ 0, and 

one can keep in the expansion of exp(— pvll) only the 
terms up to the second order. Then, Eq. (H) yields: 



167r2r2 



(6) 



where St ^ 1 — T/Tc- This is the anisotropic linearized 
GL equation of Ref. [|: 



{n\,vk) . (7) 



The anisotropy parameter in this domain follows: 



Hc2.,a 



■2 

aa 
:2 



(8) 



It should be stressed that at Tc, the anisotropy of the 
London penetration depth 7a (Tc) = \c/\a is determined 
by the same tensor (^Q'^ViVk), i.e., 7a(7c) = Jh{Tc), see 
Refs. |l| . 



or 
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The formal method. To treat arbitrary T's, we follow 
Helfand-Werthamer's routi ne ||l3| (for Eilcnberger based 
procedure see, e.g., Refs. 0, |l4[): introduce = VxUvy 
and = ± illj^ for the field along z, write vU = 
(w+II" + w~n+) /2, and use known properties of expo- 
nential operators to evaluate exp(— pvll). There are a 
few ways to proceed with actual calculation. E.g., by 
writing uj~^ = 2 dpe~'^P'^, one can sum up over ui 
in Eq. (^ and evaluate numerically the remaining inte- 
grals over p. We found it a faster procedure to follow 
Rieck and Scharnberg ||^: first to express integrals over 
p in terms of Repeated Integrals of the Error Function 
i"erfc(a;) and then numerically evaluate convergent 
sums over to. We then obtain: 




(9) 



n+m 



{v-y 



in+m 



m — n+1 



V 



erfc [x) , 



where = i (pQ / An H 11^, x = (ll}/v±)^/2(I)o/t:H and 



7;2 -I- 7)2 



To satisiy the self-consistency equation we expand ^' in 
a complete set of functions constructed by Eilenberger 
to represent various vortex lattice solutions: 



(10) 



3=0 



In fact, we do not need the explicit form of (they 
can be found in Ref. |l^); suffices it to mention that 
^0 has the known Abrikosov form, whereas the rest are 
found by applying a+ to ^fQ: a^^n — Vn+ l^'„+i and 
a^^„ = y/n'i'n^i- After some algebra, we arrive at the 
following equation: 



E E E (^^^"™> 



a;>0 j,m n—0 



\ uj>0 / j=0 



As argued in Ref. |^, the integers j and n — m should 
be taken even. When the coefficients in front of each 'i'^ 
are set zero, the homogeneous system of linear equations 
for Ck is obtained. The determinant of this system has 
to be zero, which gives an equation for H; the highest 
root is the upper critical field. 

MgB2. The band structure calculations show that the 
Fermi surface of this material consists of sheets com- 
ing from two TT-bands and two quasi-two-dimensional a- 
bands |[ |l^ . It has been shown by solving the Eliash- 
berg equations [|[ |l^ that the gap on the four Fermi 
surface sheets has two sharp maxima: Ai 1.7 meV at 
the two TT-bands and A2 ~ 7meV at the two tr-bands. 
Within each of these groups, the spread of the gap val- 
ues is small, and the gaps can be considered as constants, 
the ratio of which is nearly T independent. Experimen- 
tal estimates of the mean-free path I show that in most 
of available samples I ^ ^(0), i.e., the material can be 
considered as a clean superconductor. 

The macroscopic anisotropics such as those of 77c2, of 
the penetration depth A, of the magnetization, etc., are 
determined within a weak-coupling theory by the Fermi 
surface averages [see, e.g., Eqs. (0) and (|)]. Therefore, 
for qualitative study of these anisotropics, fine details of 
the Fermi surfaces are unlikely to be relevant. We then 
model the Fermi surface of this material as two separate 
sheets Fi and F2 with average characteristics taken from 
actual band structure results. Thus, we consider a model 
material with the gap anisotropy given by 

n(v) = J7i,2 , V e Fi,2 . (12) 

Denoting the densities of states on the two parts as Ni^2: 
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we obtain for the general averaging: 



Ni{Xi)+N2{X2) 
7V(0) 



Z.1 (Xi) + U2 {X2) , (13) 



where A^(0) is the total density of states and we intro- 
duced normalized densities of state 1^1,2 for brevity. We 
have then instead of Eq. ph : 



fi^j/i + iljf'z = 1 , t^i + 1^2 = 1 . 



(14) 



We also assume that the two parts of the Fermi surface 
have the symmetries of the total, e.g., (v)^ = where 
the average is performed only over the first Fermi sheet. 

According to Refs. ||, ^ the relative densities of states 
vi and 1^2 of our model are « 0.56 and 0.44. The ratio 
A2/A1 ~ ^2/^1 ~ 4, if one takes the averages of 6.8 
and 1.7 meV for the two groups of distributed gaps as 
calculated in Ref. ^. Then, the normalization (^4|) yields 
r^i — 0.36 and il,2 — 1.45. One should have in mind that 
the data on the gaps ratio vary from less than 3 (see, e.g.. 



and references therein) to more then 4 |18 



In the following we use averages over separate Fermi 
sheets calculated in Ref. g (vl)^ = 33.2, (i;^)^ = 42.2, 
(vl)^ = 23, and (v^)^ = 0.5 x lO^^cm^/s^. It is worth 
noting that the ratio (u^) / (u^) is 0.79 for Fi, and 46 for 
F2, whereas for the whole Fermi surface it is 1.2 . 

We further model the pieces Fi_2 by spheroids, in other 
words, we consider quadratic energy spectra £'^(fc) — 
(kl + kl)/i2ma,f.) + k^J{2m^,^), ^ = 1,2. Within 
this simple model, evaluation of averages is tedious yet 
straightforward; to get the averages given above we need 
TTia.i/rac.i = 1.3 and ma,2/'mc,2 = 0.029. 

The expansion (^0|) converges fast at high T's, but we 
need terms up to 'I'2o to stabilize the low temperature val- 
ues of Hc2.a- The temperature dependence of Hc2 for two 
principal directions so calculated is shown in Fig. |l|. The 
curve Hc2,c{T) is similar to that of Helfand-Werthamer 
p3[ , the result confirmed by the data of Ref. [l9[ On 
the other hand, Hc2,a has a slight positive curvature at 
high temperatures. As a result, the ratio Hc2,a/ Hc2,c is 
temperature dependent, as shown by the upper curve of 
Fig. p 

It is of interest to note that within a few %, solutions 
vj* of the self-consistency Eq. (|ll|) can be approximated 
by solutions of —^ff.HiHk'^ — 5", with properly chosen 
parameter ^h{T) = ^a/S,c- The reasons for this will be 
discussed elsewhere Physically, this means that the 
angular dependence of Hc2 at any T should be close to: 



c2[ 



,T) = 



c2,ab 



V7l(7> 



(15) 



where is the angle between the applied field and the c 
axis. This angular dependence has, in fact, been recently 
reported pO|. 




FIG. 1: Temperature dependence of the upper critical field 
for two principal directions. Note: each field is normalized on 
its own slope at Tc. Tc{AHc2/AT)t^. 




FIG. 2: Anisotropy ratio 7h = Hc2.ab/ Hc2,c versus T /Tc cal- 
culated with parameters for MgB2 given in the text. Dashed 
line is 7a (r) — Xc/Xab calculated as in Ref. II2I 



The drop of the Hc2 anisotropy with increasing T has 
been recorded by Angst et al. in measurements on single 
crystals of MgB2 ||2^ . Bud'ko and Canfield used a robust 
method of extracting the anisotropy of Hc2 [@ in the 
whole temperature range using the T dependence of the 
magnetization of random powders |23 . Recent specific 
heat measurements of Lyard et al. 



I and magnetiza- 
tion data of Welp et al. \ 24 on single crystals produced 
similar results. All these data show qualitatively similar 
behavior to that of the upper curve of Fig. ^. 

Physically, the large anisotropy of Hc2 at low temper- 
atures {~ 6 in our calculation) is related to the large gap 
value at the Fermi sheet which is nearly two-dimensional. 
With increasing T, the thermal mixing with the small- 
gap states on the three-dimensional Fermi sheet sup- 
presses the anisotropy down to 2.6 at Tc- 

It is worth to compare here the Hc2 anisotropy with 
the anisotropy of the penetration depth j\(T) = Xc/Xa 
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calculated within the same model of MgB? Hi 211. The re- 



[3] I. N. Askerzad e. A. Gencer, N. Guclii, and A. Kilig, cond 



suit is shown by the dashed line in Fig. 2. At Tr-. the 



mat/0201112 



anisotropy coincides with 7;^, as they should because 
the GL theory contains only one "mass tensor" which de- 
termines anisotropics of both and A^. The expression 
for the anisotropy ratio clearly amplifies the contribu- 
tion of Fermi surface parts with the large gap. However, 
in a clean material at T = 0, the superfluid of Cooper 
pairs is Galilean invariant; this implies that all charged 
particles participate in the superflow, their energy spec- 
trum notwithstanding. This is reflected in the result for 
the A-anisotropy at T = 0, 



7^(0) 



(16) 



neither the gap nor its anisotropy enter this expres- 
sion. As mentioned, for MgB2, the ratio (|l^ ) is 1.2, i.e., 
7a(0) = 1.1 in a striking difference with the large zero-T 
anisotropy of Hc2- 

Thus, the question "what is the anisotropy parameter 
of MgB2?" does not have a unique answer. To pose the 
question properly one should specify the quantity of in- 
terest. If this is Hc2, the answer is given by the upper 
curve of Fig. 2 for a clean material; if this is A, see the 
dashed line. If this is the magnetization in intermediate 
fields, M (X ((/)o/A^) \n{Hc2/ H), the main contribution to 
anisotropy comes from A, however, the Hc2 anisotropy 
contributes as well (being smoothed by the logarithm). 
The last situation should be taken into account while ex- 
tracting anisotropy from the torque data in tilted fields 
p5| , p6[ , the point stressed by Angst pTf . 

Given the simplifying assumptions about the Fermi 
surfaces we have made and our results for the Hc2 
anisotropy which reproduce qualitatively well the mea- 
sured anisotropy, we believe that the anisotropic prop- 
erties of Hc2 and A we have described, are generic for 
materials with anisotropic gaps and Fermi surfaces. 
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